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Abstract 

The resummed Drell-Yan cross section in the double-logarithmic approximation suffers 
from infrared renormalons. Their presence was interpreted as an indication for non- 
perturbative corrections of order Aqcd/(Q(1 — z)). We find that, once soft gluon emis- 
sion is accurately taken into account, the leading renormalon divergence is cancelled by 
higher-order perturbative contributions in the exponent of the resummed cross section. 
From this evidence, 'higher twist' corrections to the hard cross section in Drell-Yan pro- 
duction should intervene only at order Aq CD /((Q 2 (1 — z) 2 ) in the entire perturbative 
domain Q(l — z) > Aqcd- We compare this result with hadronic event shape variables, 
comment on the potential universality of non-perturbative corrections to resummed cross 
sections, and on possible implications for phenomenology. 
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1 Introduction 



The factorization theorems of QCD [|]] allow a rigorous treatment of 'hard processes', 
since the non-perturbative dynamics can be isolated in a few universal distribution and 
fragmentation functions. The classical case for this approach is the Drell-Yan (DY) 
process h a + hb — > /i + ^ + X, when a muon pair (or, alternatively, a massive vector 
boson) with invariant mass Q 2 is produced in the collision of two hadrons, h a , hb, with 
invariant mass s. The cross section is then given byQ PJ 

a W(s,Q ) a = ——^ — (1.1) 



dQ 2 9 Q2 S 
1 / k\ 

W{s,Q 2 ) = Y,J ^^h/ ha {^Q 2 )f 3 /h b (x 3 ,Q 2 )^{z ias (Q)) + Of (^) ) , (1-2) 

where fi/h a { x h Q 2 ) (fj/h b { x ji Q 2 )) is the distribution function for a parton z in h a (parton 
j in hb) and z = Q 2 /(xiXjs) is the ratio of invariant masses of the produced muon pair 
and the colliding hard partons. The 'hard cross section' u>ij(z,a s (Q)) can be expanded 
as a power series in the strong coupling a s (Q). The separation of the Drell-Yan cross 
section into distribution functions and a hard cross section is not unique. We will take 
the DIS scheme for the distribution functions, so that for quarks and antiquarks, which 
is the only case we will be interested in, the distribution functions are given by the 
corresponding deep inelastic scattering (DIS) cross section. In the following we neglect 
the sum over parton species and consider only i — q, j — q. 

The DY process is one of the simplest, where one can realize, experimentally, a short- 
distance dominated dynamics involving two large but possibly disparate scales: When 
1-z C 1, the phase space for gluon emission is restricted, so that the actual scale of the 
emission process is Q(l — z) rather than Q. Although soft divergences cancel, large finite 
corrections in the form of '+- distributions' a s (Q) n [(ln 2n " 1 (l — z))/(l — z)} + are left over 
and it is necessary to resum them to all orders. As long as Q(l — z) > Aqcd ; the relevant 
contributions from soft and collinear partons are still amenable to perturbative analysis. 
The resummation (to single-logarithmic accuracy) has been completed in @, 0| . 

The theoretical discussion is simplified in terms of moments. Introducing the DY 
scaling variable r = Q 2 /s, eq. ( |1.2| ) can be written as 



i 

W(N,Q 2 ) = J dTT N - l W(T,Q 2 ) = f q/K {N,Q 2 )f- q/hb {N,Q 2 )^{N,a s {Q)), (1.3) 
o 

where u qq -(N, a s (Q)) (f q / ha (N, Q 2 )) denotes moments of the hard cross section (distribu- 
tion function) with respect to z (xi). For large N, the moments probe small 1 — z with 
the correspondence rule [(ln n (l — z))/(l — z)] + < — > ln n+1 N. The large higher order 
corrections a s (Q) n \n h N (k < 2n) can be resummed systematically by 0, 0] 

1 For simplicity wc quote the Born cross section for production through an intermediate photon only, 
and neglect the contribution from Z . 
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u g q(N,a s (Q)) = H(a s (Q))exp 



JV-1 



dz 



1 - z 



dk 2 t 



+ B(a s (VT^Q)) 



+ o 



kiN 



(1.4) 



provided the functions A(a 8 ) and B(a s ) (to be detailed later) are known to some appro- 
priate order. 

In this paper we address the question of possible non-perturbative contributions to 
the resummed cross section in eq. ( |1.4| ), suppressed by powers of the large momentum 
Q ('power corrections'). We are motivated by the observation that the integrals in this 
equation include the regions of very small momenta k t , which render the resummation 
of In N sensitive to the infrared (IR) behaviour of the strong coupling. The resummed 
expression depends on the prescription how to deal with the IR Landau pole in the 
(perturbative) strong coupling or, in other words, on the value of an IR cutoff. It has 
been found j?J |5| that IR effects to the resummed formulas should intervene already at 
the order Aq CD /(Q(l — z)), which suggests that to this accuracy the factorized hard 
cross section has to be complemented by (exponentiated) non-perturbative corrections. 

If confirmed, this statement is of a considerable theoretical and phenomenological 
interest. For phenomenology, it would imply that non-perturbative corrections decrease 
so slowly that they are numerically important for data analysis even at the largest avail- 
able energies. For a theoretician, the problem is especially interesting because without 
guidance of the operator product expansion there is no immediate classification of power 
corrections and one must rely on different methods. 

The approach applied in [|], [5] and in this paper utilizes the fact that the insufficiency 
of perturbation theory to account for the region of momenta of order Aqcd is indicated 
by the perturbative expansion itself, through the appearence of another type of large 
perturbative corrections in large orders. These corrections, proportional to a s (Q) n n\ 
and known as infrared (IR) renormalons ||, cause the divergence of the perturbative 
expansion. The numerical value attached to such a divergent series depends on the sum- 
mation prescription. Since, according to current understanding, QCD is the theory of 
strong interactions, this prescription dependence must be cancelled by all those contri- 
butions that elude a perturbative treatment. By this argument one can determine the 
order in Aqqb/Q at which non-perturbative corrections must enter, or, in other words, 
the level of theoretical accuracy beyond which perturbative QCD must fail. For pro- 
cesses that allow an operator product expansion like in DIS, the renormalon approach is 
complementary and less powerful than the operator product expansion. In the case of 
resummed cross sections the consideration of IR renormalons or, equivalently, IR cutoff 
dependence in perturbation theory can provide genuinely new information on the na- 
ture of non-perturbative corrections. The results of (|, |[ suggest that non-perturbative 
corrections are expected to affect the resummed DY cross section at the level of 1/Q. 
Corrections of this size were also detected for hadronic event shape variables f7|, || . Since 
the evolution equations that govern soft gluon emission are universal, it was proposed 
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that non-perturbative corrections of order 1/Q can be described by a single parameter 
common to all resummed cross sections |s|, [Top- 

In this paper we reanalyze this problem. Our main result can be summarized by the 
statement that the structure of IR renormalons and power corrections depends crucially 
on soft gluon emission at large angles so that the collinear approximation does not 
apply The phase space is then much more complicated and pro cess- dependent. To 
obtain information on soft gluon emission to power-like accuracy, the functions A and 
B in eq. ( |1.4|) need to be kept to all orders. In particular, the leading IR renormalon 
divergence (implying 1 / Q-corrections) of the resummed DY cross section found in ||. |3| 
appears as an artefact of using a finite-order approximation for these functions. The 
apparent l/Q-ambiguities in the 'standard' factorized formula eq. flLJ ) are cancelled a 
zero in higher order perturbative corrections and the remaining power corrections appear 
to be of order N 2 Aq CB /Q 2 . This cancellation is in effect similar (though the physics is 
different) to the cancellation of the leading 1/mg ambiguities between the phase-space 



and the series of radiative corrections to the decay widths of heavy particles [11], |T^, [13]. 
Our results are obtained in a certain approximation to large-order behaviour. To go 
beyond this approximation, two-gluon emssion would have to be analyzed to power-like 
accuracy. The fact that even for one-gluon emission the correct IR renormalon structure 
follows only from all-order expressions for A and B makes this a vastly complicated 
problem. 

Since the language connected with large-order behaviour appears rather formal it 
is helpful to recall the following correspondence which reveals more physical insight: 
For processes to which the operator product expansion can be applied, the power cor- 
rections deduced from renormalons can be identified with contributions of operators of 
higher dimension (twist). The above prescription dependence translates into the ambi- 
guity in the choice of factorization scale below which fluctuations are called 'soft' and 
therefore should be included into the definition of the matrix elements of higher-twist 
operators. For perturbative contributions this factorization scale acts as an IR cutoff. 
Therefore, more generally, quite in analogy to the large perturbative corrections of type 
a s (Q) n [(\n 2n ^ 1 (1 — z))/(l — z)] + discussed previously, the large factorials a s (Q) n n\ can 
also be understood [12 as remnants from the Kinoshita-Lee-Nauenberg infrared cancel- 
lations, provided an IR regularization with an explicit mass scale such as a finite gluon 
massf] A is used. Instead to factorials in large orders, one can then restrict attention 
to cutoff dependences like X/Q, X 2 /Q 2 \nX 2 in lowest order radiative corrections, which 
is physically more intuitive. For IR safe quantities, explicit IR divergences vanish by 
construction, while the scaling of power corrections can be determined by setting the 
IR cutoff to Aqcd- In this language, we find that potential contributions to Drell-Yan 
production of order X/Q from the restriction on the phase space are exactly cancelled 
by the contributions of the same order to the matrix elements.^ 

2 The precise correspondence of fl2|| is true in a certain approximation, where only diagrams without 
gluon self-coupling are considered, so that no difficulties with gauge-invariance occur. A generaliza- 
tion, although physically compelling, is unknown, reflecting the present ignorance how to deal with 
renormalons on a diagrammatic level, when inclusion of gluon self-couplings is necessary. 

3 More precisely, we find that the matrix elements can not be expanded in powers of X/Q near the 
phase space boundaries, and should be treated exactly. 
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The paper is organized as follows: In Sect. 2 we analyze the exponent in eq. ( |1.4j ). We 
reproduce the results of [|, || in the approximation adopted in these works and proceed 
to show that low-order truncations of the functions A and B are insufficient to draw 
definite conclusions about non-perturbative corrections to the exponent. We outline 
three possible scenarios for the complete result. To distinguish the correct one, an all- 
order calculation beyond the soft and collinear approximation is required. We introduce 
and justify the corresponding (still approximate) all-order calculation in Sect. 2.3. Those 
readers who prefer to think in terms of IR cutoff dependence rather than large-order 
behaviour might ignore subsection 2.3. 

In Sect. 3 we calculate the partonic DY and DIS cross sections with finite gluon 
mass A to first order in the strong coupling and extract the hard scattering function 
cu q q(z,a s (Q)). We find that all potential corrections of order NX/Q cancel, so that the 
one-loop perturbative result is protected from IR contributions to this accuracy. As 
mentioned previously this calculation has a dual interpretation in terms of either cutoff 
dependence or large-order behaviour. At this stage we can conclude that at present 
there is no evidence for nonperturbative corrections of order NAqqd /Q to the DY cross 
section. To trace the origin of this cancellation, we repeat the calculation in the soft limit 
in Sect. 4, but with a cutoff on the energy and transverse momentum of emitted gluons. 
This allows us to pin down the phase space approximation, in which the DLA looses 
terms of order NAqqb/Q arid clarifies the relevance of soft emission at large angles. 

In Sect. 5 we follow the approach of [|14j to derive the soft factorization for the DY 
cross section in terms of Wilson lines. We calculate the corresponding Wilson line with 
an arbitrary number of fermion loop insertions in the gluon line and check that the result 
satisfies the correct renormalization group equation (evolution equation). We find that 
all relevant anomalous dimensions (in the MS scheme) are entire functions in the Borel 
plane and the IR renormalons enter only through the initial (boundary) condition for 
the evolution. The apparent l/Q-corrections to the standard result eq. ( |L4p are due to 
an unfortunate choice of particular solution of the evolution equation for the Wilson line 
which in turn necessitates a more singular homogeneous solution than actually required. 
As a result, the apparent IR renormalon that indicates l/Q-corrections in the Borel 
transform of the exponent in eq. (|1.4fl is cancelled by a 'hidden' zero in the sum of 
Borel transforms of the (redefined) anomalous dimensions, which becomes manifest only 
if these are calculated to all orders. We discuss a possibility to reformulate the standard 
factorization formulas for DY cross sections in such a way that the cancellation of leading 
renormalons is explicit. 

In Sect. 6 we summarize, compare our results with those for thrust averages in e + e~ 
annihilation and comment on potential implications for processes other than Drell-Yan. 
Appendix A contains a rederivation of the cusp anomalous dimension of the Wilson line, 
and in Appendix B we compare typical phase-space integrals for the DY cross section 
and the thrust distribution. 
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2 Anatomy of the exponent 



In this section we discuss the exponent in the resummed hard cross section. We rewrite 
eq. (p~SD as 



u q g(N,a s (Q)) = H{a s (Q)) exp[E(N,a s (Q))} + R(N,a s (Q)) , (2.1) 
where R(N, a s (Q)) vanishes as N — > oo and the exponent is given by 

2 (l-z) 2 

dfc? 1 




E(N,a s (Q)) = Jdz — 1 2 J -± A(a s (k t )) + B(a s (VT^Q)) + C(a s ((l- z)Q)) ^ . 

(2.2) 

The exponent expanded in a s (Q) contains terms a s (Q) n \n h N with k < n + 1. To sum 
all logarithms with k > n+l — m, the expansion coefficients of A(a s ) = yj m=0 a n a™ +1 U P 
to order m + 1 and of B and C up to order m are required, as well as the /3-function that 
controls evolution of the coupling to order a™ +1 . A is a process-independent function, 
which can be identified with the eikonal (or cusp) anomalous dimension. Up to second 
order it is given by [ |T5p 



Cp Up 
it 2ir z 




7T 2 \ 5 Af 

— - -N f 
6 9 s 



(2.3) 



B and C are process-specific. B involves the underlying DIS process, while C comes 
only from the DY process. Their first order expressions are 



c = 0. (2.4) 



3C F 

Note that compared to the conventional form quoted in eq. (|1.4|) we have reintroduced 
a third function C with expansion in the coupling normalized at (1 — z)Q as it was 
originally present in the analysis of . With the coefficients given above, C is irrelevant 
for summing logarithms to next-to-leading accuracy (k > n) and can be dropped for 



this purpose Jl6|, |17]]. The reason for keeping this function in the present context will 



become clear in Sect. 5. The function C can be dispensed of by a redefinition of A and 



B ||17|| . However, the interpretation of A as universal cusp anomalous dimension is then 
lost beyond order a 2 s . 

We are now interested in the infrared sensitivity of the exponent. Perturbatively, 
it is most directly visible through the presence of the Landau pole of the perturbative 
running coupling inside the integration regions in eq. ( |2.2| ). The Landau pole has two 
(related) consequences. It leads to IR renormalons in the re-expansion of E in a s (Q) and 
renders the integrals dependent on the treatment of the pole. The ambiguities that arise 
in this way quantify our ignorance on the correct infrared behaviour, which does not 
show a Landau pole. Since extracting the ambiguities from the Landau pole essentially 
projects on the infrared regions of integrals, the same conclusions can be obtained from 
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eliminating these regions by an explicit IR cutoff and by studying the dependence on 
the cutoff. In Sect. 2.1 we reproduce the result on power corrections in [|J obtained in 
the double-logarithmic approximation. The general case is developed in Sect. 2.2. We 
introduce and justify an approximation to calculate large-order behaviour in Sect 2.3. 



2.1 Ambiguities to double- logarithmic accuracy 

The double-logarithmic approximation (DLA) corresponds to exponentiation of the 
a s (Q) In 2 iV-term in uj q q{N 1 a s {Q)). In eq. Q2.2p, only a has to be kept.Q We remove 
the contribution from gluons with energy less than /i > Aqcd- Since Q{1 — z)/2 can be 
interpreted as the energy of the emitted gluon (for small 1 — z), we obtain 



E DLA (N,a s (Q),v) = J dz i _ - 2a J 



2(1-*) % 



/i-independent + 2a a s (Q) {N - 1)^ In ^ + O (j^J . (2.5) 



Notice the linear dependence on the cutoff and the large coefficient proportional to N. 
In [|j] the ambiguity of the exponent due to the Landau pole has been evaluated in the 
approximation of one-loop running for the coupling, 

,,s = a s (Q) = 1 (2 fi s 

as{ ] l-/3 a s (g)ln(fc 2 /Q 2 ) (-/3 ) M^/A 2 CD ) ' l " J 

where (3q = — 1/(4tt)(11 — 2Nf/3) is the lowest-order coefficient of the /^-function. Fol- 
lowing the approximation of we keep only the term aoa s (k t ) up to terms that give 
ambiguities of order l/Q 2 . The fc r integral produces a cut starting at z — 1 — Aq CD /Q. 
Evaluating the ^-integral above or below the cut, we find the difference (neglecting terms 
of order (A QC d/Q) 2 ) 



1 



y N-l 



5E ULA (N,a s (Q)) = —Jdz 

o 

= 2a (iV - 



— 2a f-TT ) disc 

z V Po 



1 \ A, 



QCD 



PoJ Q 



lnln(l 



Q- 



' NA 



A 2 



QCD V 



Q J 



(2.7) 



The cutoff dependence of eq. Q2.5| ) and the ambiguity in eq. ( p.7| ) agree, if we choose 
fi ~ Aqcd, since a s (Q) ln(AQ CD /Q 2 ) = (— 1/fio). It can be shown |l8j that the ambiguity 
does not change qualitatively, if the coefficients a\ and bo are included as required for 
the summation of next-to- leading logarithms a s (Q) n ln n N. We therefore conclude that 
with the truncated series for A and B sufficient to sum next-to-leading logarithms, the 



4 We use the term DLA, if only ao is kept, irrespective of whether a s that multiplies ao is frozen or 
running. 
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exponent in the form of eq. ( |2.2|) is prescription-dependent by terms of order NAq C b/Q, 
if the expression for the running coupling is not expanded and truncated. At this point, 
however, it is not clear whether these terms are artefacts of the approximations made or 
whether they should be interpreted as an indicator for 'higher twist' effects of this order 
as suggested in ||, |5|]. To illustrate this point, we note that to the accuracy of summing 
next-to-leading logarithms one may replace || 

z N ^ - 1 - -G (l - - z) (2.8) 

It is easy to check that after this substitution the corresponding exponent E has no ambi- 
guities at all, unless A" > Q/Aq C b- Such large moments are excluded from consideration, 
since there is no perturbative treatment for them.0 

2.2 General case 

In higher orders (contributing only to sub-dominant logarithms a s (Q) n \n k N, k < n), 
the replacement in eq. (|2.8|) requires a redefinition of the functions B and C. This sug- 
gests that for the purpose of identifying power corrections to the exponent, higher order 
coefficients of A, B and C are relevant, while they are not for a systematic summation 
of logarithms up to a certain accuracy. The discussion of ambiguities including higher 
order coefficients is simplified in terms of Borel transforms. We write the exponent as 

oo 

E(N,a s {Q)) = (^-j^j J d«er«/(-/WQ)) B[E](N,u), (2.9) 

o 

where for any function 



f(a s (Q)) = ]T f n a s (Q) n+l B[f](u) = £ /„(-#,)-»_ . (2.10) 

n=0 n=0 n ' 

Note the expansion parameter is a s (Q). If / has a term of order a s (Q) , it is treated 
separately. The ambiguities of E due to the pole in the running coupling are now 
generated by IR renormalon poles of B[E] in the u- integral. (These poles correspond to 
factorial divergence of the corresponding series expansion in a s (Q).) A pole at u = m 
leads to an ambiguity of order (Aq CD /(5) 2m . The Borel representation allows us to write 
B[E] in a factorized form. Let us continue to assume that f3(a s ) = Poct^, so that eq. 
is exact. Then 



00 ( 1 \n n — 1 

, —L\nx 2 u n = x- 2u B[a s (xQ) n ](x,u) = 
' n\ n! 



B[a s (xQ)}{x,u) = 22 ^d~ lnx un = x~ B[a s (xQ) n ](x,u) = -^-x~ 2u , (2.11) 

n=0 



5 It can be seen that the ambiguities that arise for N > Q/Aqcd in this case are not related to IR 
renormalons (poles in the Borel transform, see below), but to the divergence of the Borel integral at 00, 
when one of the kinematic variables, here Q/N, becomes of order Aqcd, see Jis| . 
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so, for any function f(a s (xQ)), 

B[f(a s (xQ))](x,u) = x~ 2u B[f{a s {Q)](u) . (2.12) 

The k t - and z-integrals that appear in E can then be done and we obtain 



B[E](N,u) = 2A 1 (N,u)B[A](u) + A 2 (N,u)B[B](u) + A 3 (N,u)B[C](u) 



(2.13) 



where all dependence on N is factored into the functions 



A '( jv '"»=Ki +r < jv ) 
A2(AVtl)= i + r < w > r (-*> 



T(-u) T(-2u) 



u T(N-u) ' 



T(N-u) T(N-2u)_ 

A 3 (N,u) = A 2 (N,2u) , 



(2.14) 
(2.15) 



which depend only on the form of the exponent, but not on the details of the functions 
A, B, C. Note that A 3 = A 2 - uAi, so that 



B[E]{N,u) = A 1 {N,u){2B[A]{u) -uB[C](u)} + A 2 (N,u) [B[B)(u) + B[C](u)] . (2.16) 



This is the observation that C is in fact redundant |T7|] and can be absorbed into a 
redefinition of A and B. If, as in Sect. 2.1, we keep only the coefficient a$ of A, we 
obtain close to u = 1/2, 

B[E DLA ](N,u) u ^ /2 2a (N-l)—^—. (2.17) 

1 — 2u 



The ambiguity of E derived from eq. (|2.9|) (evaluating the integral above and below the 



pole) coincides with eq. (|2.7j ). Without truncation of A and C, we obtain 



B[E](N,u) u ~"= 2 2(N - 1) ' 



l-2u 



25 [A] ( ^ 



(2.18) 



assuming only that -B[A] and B[C] are non-singular at m = 1/2. Notice that the term 
A2B[B] and the first term in square brackets of Ai which originate from the DIS subpro- 
cess do not lead to a singularity at u — 1/2. The leading pole is at u — 1, indicating power 
corrections of order NAq CD /Q 2 , in agreement with the operator product expansion for 
deep-inelastic scattering, which constrains corrections to the leading-twist approximation 
to enter at order l/Q 2 . 
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Hence, a conclusion regarding the presence of a pole at u = 1/2 can not be obtained 
from finite-order approximations^ of A etc., since the residue involves the Borel trans- 
forms at the point u — 1/2. In other words: The presence of nonperturbative corrections 
of order NAqcb/Q can only be discerned, if the anomalous dimension functions in the 
exponent are themselves defined to power-accuracy. Since the whole point of the resum- 
mation formula eq. ( |1.4|) is to obtain certain large contributions to u q q(N,a s (Q)) to all 
orders by finite order calculations of A etc., there is no compelling motivation to restrict 
attention to E alone instead of the complete u q q(N, a s (Q)) for the purpose of identifying 
power corrections, because finite order approximations to the exponent are insufficient. 
But since the exponent captures all leading singular contributions ln fc N (but not singu- 
lar terms 1/N\n k N), one would still expect the exponent to indicate also the leading 
power corrections (NAqcd/Q) 11 (but not l/N(NAQcr>/Q) k )- Regarding the presence of 
iVAQCD/Q-corrections, the following three scenarios can be envisaged: 

A. A renormalon ambiguity of order NAqqd/Q, i.e. a pole at u — 1/2 is present. 
This situation is realized if IB [A] (u) — uB[C](u) is non-zero at u = 1/2 and 
provided no cancellations with the remainder R(N,a s (Q)) occur, see C below. 
Then a nonperturbative correction of order NAqqd/Q is required, and it has to 
exponentiate. This scenario is suggested by @. ||] . 

B. The pole at u = 1/2 cancels because 

2B[A](^j-^B[C](^j=0, (2.19) 

and no ambiguity of order NAqqd/Q remains. Then no nonperturbative correc- 
tions are required to this accuracy. Note that the Borel transforms of A, B and C 
could themselves have poles. Even in this case B can not participate in the can- 
cellation (if it occurs) at u = 1/2, because A 2 (A r , 1/2) ~ N 1 / 2 , while the residue of 
the pole at u = 1/2 is of order N. Note also that after elimination of the redundant 
C to arrive at the conventional eq. (|1.4|) , the previous equation implies that the 
Borel transform of A as it appears in eq. ( |1.4| ) has a zero at u — 1/2. We will verify 
(in a certain approximation) that this zero is indeed present. 

C. A further possibility is that the exponent does not incorporate all leading (in N and 
1/Q) power corrections and the pole is cancelled with the remainder R(N,a s (Q)) 
in eq. ( [2.1|) . This can happen if the series for the remainder looks like (for large n) 

1 n 1 

R(N, a s (Q)) = ]T a s (Q) n+1 (-2(3o) n nl - £ - ln fc N 2 . (2.20) 
n iV fc=o K - 

6 This situation must be distinguished from the more familiar statement that the overall normalization 
of renormalon divergence is not obtained correctly from diagrams with a single chain of loops as discussed 
in Sect. 2.3. In the present case, even for a single chain, the overall normalization is incorrect with any 
finite-order approximation to A. 
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Figure 1: af-contribution to the partonic Drell-Yan cross section. 7* represents a 
photon with invariant mass Q 2 that splits into a lepton pair. 



To any finite order in a s (Q) the remainder vanishes as N 
transform 



00, but its Borel 



B[R](N, u) 



N 



4u-l 



1 - 2u 



(2.21) 



has a pole at u — 1/2 with residue proportional to N. A cancellation of this type 
would again imply that no nonperturbative corrections of order 1/Q are required, 
but that the soft factorization techniques break down to power-like accuracy. 

To single out the correct scenario and clarify the size of expected nonperturbative 
corrections, u q g(N, a s (Q)) or at least the functions A, C would have to be calculated to 
all orders. Obviously, this requires a certain approximation. In the following sections, 
we calculate a subset of higher order corrections and argue that this subset provides 
some generic insight. Before introducing this subset in the next subsection, we comment 
on the approximation of neglecting higher orders in the /3-function in the derivation 
of eq. (|2.13|) . When the next coefficient, f3i, is retained a simple factorized form as 
in eq. ( |2.13p is not immediate. However, provided the /^-function itself does not have 
renormalons (as believed in the MS-scheme), inclusion of higher orders are expected only 
to turn the poles in u into branch point singularities. The essence of the arguments above 
remains unaffected. 



2.3 Higher order approximation 

In calculating higher-order contributions to the DY hard cross section, we resort to an 
approximation that - despite its largely heuristic character - has proven useful in iden- 
tifying IR renormalons and power corrections semi-quantitatively in other applications 
[p0~[ |2~T| . The approximation consists in selecting the subset of higher- loop diagrams gen- 
erated by insertion of an arbitrary number of fermion loops in the gluon line in diagrams 
contributing to the lowest order radiative correction, see Fig. [I] for an example. The 
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corresponding series in a s (Q) grows indeed as a s (Q) n+1 n\ and its Borel transform has 
the expected IR renormalon poles. The approximation ('single- chain approximation') 
has two immediate deficiencies which we discuss first. One may keep in mind that even- 
tually the best justification for the single-chain approximation derives from its relation 
to IR cutoff dependence, which is discussed towards the end of this subsection. 

The first observation is that each fermion loop is proportional to the abelian part 
of /3q, i.e. Nf/(6ir). The large-order behaviour from infrared regions is sign-alternating 
and does not lead to any ambiguities. Evidently, this happens because with fermion 
loops alone, there is no IR Landau pole in the running coupling. The usual way to deal 
with this deficiency is to restore the non-abelian /3 by hand.[] This is justified by the 
expectation that the factorial behaviour is related to evolution of the coupling so that all 
other uncalculated diagrams would combine to reproduce this ad-hoc manipulation, up 
to an overall normalization, which is of minor interest for the parametric size of power 
corrections. The single-chain approximation then amounts to integrating the lowest order 
corrections with the complete (but still approximate) gluon propagator a s {k)/k 2 . This 
deficiency is not specific to the DY process, but common to all previous applications 

mm- 

A more specific draw-back of the single-chain approximation is that within the se- 
lected set of diagrams, exponentiation of soft and collinear logarithms IniV does not 
occur. The calculated diagrams are those with the largest number of factors Nf (the 
number of massless fermions) and give at most a s (Q) n ln n+1 N (This remains obviously 
true after restoring the non-abelian /3q.). The dominant terms in the large- N limit, 
a s (Q) n ln fc N (n+ 1 < k < 2n, n > 1), come from diagrams with two or more gluon lines 
(chains). They can be partially obtained by exponentiation of the single-chain result, but 
in any case exponentiation leads outside the set of diagrams calculated exactly. Thus, if 
u q q (N,a s (Q)) denotes radiative corrections in the single-chain approximation (SCA), 
normalized to the tree-level hard cross section, 



In 



1 + J q S q C) (N,a s (Q))} = J q S q C) (N,a s (Q))+ terms beyond SCA. (2.22) 



These terms are of order ((In N)/Nf)) k relative to uj qq ^ C \N,a s (Q)) so that the large-iV 
limit (the one we are interested in) does not commute with the large- Nf limit (the one 
we calculate). 

At this point, we note that the summation of soft and collinear IniV's has two as- 
pects: The first is exponentiation as a consequence of an evolution equation, the second 
the arguments of the couplings that appear in the exponent, eq. (|2.2j). In the present 
investigation of power corrections to the DY process, as well as all previous ones [[|, |5| , 
we are mainly interested in the exponent and therefore the second aspect, which is in- 
dependent of whether the evolution equation is an equation for In [u qq (N, a s (Q))] or for 
id qq -(N, a s (Q)) (in which case exponentiation does not occur). Indeed, the single-chain 
approximation gives a non-trivial series of higher-order corrections to the anomalous di- 
mension functions A, B and C from subdominant logarithms in diagrams with fermion 



7 This is reflected already in our definition of the Borel transform, eq. (2.10), where /3o refers to its 
non-abelian value. 
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loops. The fact that this series exponentiates beyond the adopted approximation is 
secondary as long as one is interested only in the consequences of integration over the 
running coupling in the exponent. One could also say that it is more appropriate to think 
of this approximation as an approximation to ln[uj qq -(N,a s (Q))]. When interpreted in 
this way, corrections are formally of order 1/Nf with no factor of IniV and the above 
non-commutativity of the large- N and large- Nf limit is, at least superficially, absent. 

In the context of consequences of integration over the running coupling, it is more 
important that the single chain approximation is consistent with the fact that the ap- 
propriate scale in the coupling to sum large logarithms close to threshold is given by the 
transverse momentum of the emitted parton (gluon) |22| 



The real emission correction 



to the partonic DY process in this approximation is schematically given by 



W ie&1 (z,Q'< 



dk 2 

disc fc2 / -p- 



d\k\d 2 k t 



F(z;\k\,k t ;k 2 )a s (k), 



(2.23) 



(2tt) 3 2\/ k 2 + k 2 



Kinematics dictates that the discontinuity is integrated up to (k 2 



Q 2 (l - z) 



Interchanging the order of the k - and kt- integral, k is integrated up to A; 2 < Q (1— z) /4 



(for z close to 1). Following the argument of [22| the fc 2 -integral can be evaluated up to 
subleading logarithms in 1 — z by setting k 2 to the value at the upper limit of integration 
in all terms singular as k 2 — > and to zero elsewhere. Therefore 



W rea \z,Q 2 ) ~ 



d 3 k 



(2tt) 3 2|A:| 



^-F(z;\k\,k t ]0) a s (kt), 



(2.24) 



up to subleading ln(l — z)'s. The previous line gives simply the one-loop correction, 
but with a s normalized at kt ~ Q(l — z). We see that the single-chain approximation 
is fully consistent with p2|: The only approximation is that the the gluon-propagator 



that enters the Schwinger-Dyson equation of [^] is given by a s (k)/k 2 in the single-chain 
approximation. In this approximation, the argument of the coupling is determined by 
the subprocesses in which the gluon line splits into the qq pair, with k 2 being the largest 
possible invariant mass of this pair, which is allowed by the kinematics of soft emission. 

Let us emphasize that due to the approximations made to arrive at eq. ( f2.24| ), it 
can not be used immediately to deduce power-corrections. Also, there is no justification 
for integrating eq. ( p.23| ) with a s (k t ) instead of a s (k), since there is not even a formal 
limit in which such a replacement has a diagrammatic interpretation or would model a 
complete gluon propagator. 

It is useful to observe that the Borel transform of radiative corrections in the single- 
chain approximation can be calculated as [O, 



B 



UJ 



(SC) 

<?<? 



(N,u) 



sin(7ru) 



TTU 




d\ 2 



^(N,X 2 /Q 2 ) 



(2.25) 



where u q g (N,\ 2 /Q 2 )a s (Q) denotes the one-loop correction to the hard cross section, 
calculated with a gluon of mass A. The factor 5/3 in the exponent arises, because we 
renormalize fermion loops in the MS scheme. 
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Eq. ( p. 25 ) gives a transparent interpretation of IR renormalons in terms of IR cutoff 



dependence. By general properties of Mellin transforms, poles at positive u on the 
left hand side of ( [2.25 ) are related to non-analytic terms (in A 2 ) in the expansion of 



oj(fi(N, X 2 /Q 2 ) for small A 2 . The existence of the pole in question, at u = 1/2, is related 
to the presence of a/A 2 , at u = 1 to A 2 In A 2 etc.. Indeed, instead of examining the large- 
order behaviour of perturbation theory to find hints for power corrections, it appears 
more direct to explore the IR cutoff dependence in lowest order. Then, just as an IR 
divergence In A 2 in the partonic DY cross section indicates that the process is not IR 
safe but requires the introduction of IR sensitive distribution functions, the presence of 
further non-analytic terms indicates 'higher twist' nonperturbative corrections, at least 
with the adopted choice of leading-twist distribution functions. 

Eq. ( p. 25 ) is technically convenient, because to calculate the hard cross section an 



IR regulator for intermediate steps is needed even in any case. Taking finite gluon mass 
instead of the more conventional dimensional regularization, the calculation differs from 
the one-loop calculation only in that we do not neglect terms that vanish as A — > 0. 

3 Hard cross section with IR cutoff 

In this section we calculate the hard Drell-Yan cross section with finite IR regulator and 
investigate its cutoff dependence. The choice of gluon mass A as regulator is convienent 
because by eq. ( 2.25|) the cutoff dependence can be translated into a statement on IR 



renormalons. 

The calculation is a textbook problem. Since u qg -(z,a s (Q)) does not depend on 
the particular initial state, one chooses to calculate the DY cross section and structure 
functions for quark states. It is convenient to introduce the dimensionless variable £ = 
X 2 /Q 2 . For the partonic DY and quark distribution function, let us write 

W(r,Q 2 ,0 = Q\ [8(l-T) + W^(r,0 + O(a s (Q) 
f q / g (x,Q 2 ,0 = 5(l-x) + f^ q (x,0 + o(a s (Q) 2 ), (3.1) 

where Q q is the electric quark charge in units of electron charge. In the DIS scheme the 
distribution function is defined as 

M x > Q 2 ) = f = ~K + 12x 2 ^fW^ . (3.2) 

W^ v denotes the hadronic tensor for deep inelastic scattering from quarks with momen- 
tum p. According to eq. ( p..2|) , the hard cross section is given by 

u qq {z,a s (Q),0 = Q 2 g [6(l-z) + uV(z,0 + o(a s (Q) 

4 } M) = wW(z,0-2/$(*,0, (3.3) 

using equality of /,/, and f q/q . 
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3.1 Partonic Drell-Yan cross section 

The matrix element for 



g(pi) +q{p 2 ) — >g(k) +7*(g) 



(3.4) 



averaged (summed) over initial (final) colours and polarizations is given by 



\M\ 



DY 



, Cf9 2 s 



2s(Q 2 + X 2 ) t X 2 Q 2 u \ 2 Q 2 



tu 



+ - 



u u 



2 + -t 



t 2 



(3.5) 



with s = (pi +P2) 2 = Q 2 /z, t = (pi — A;) 2 , u = (j>2 — k) 2 . Integrating over the momentum 
of the photon with invariant mass Q 2 gives 



wW> Teal {z,o 



and, finally (z = 1 — z), 



d 3 k 



J 2 v^T^ 



s((pi+P2-k) 2 -Q 2 ) \M\ 2 DY , (3.6) 



7T I V Z — £ Z 



* \ z-tjz + ^/{z + tjz)'> - Atjz 

, — V z — tz In . 

\ V z - & s y z-iz- y/(z + izf - 4£z 



(3.7) 



The virtual gluon corrections are 



W^\z,0 = ^^8(l-z){(l + 2 \u 2 (-0-l:ln 2 t + lnZHl + + t 

7T I 2 fa 



3 7 

2^-4-^-* 



(3.8) 



The double-logarithmic divergences In 2 £ from soft and collinear gluons cancel between 
real and virtual corrections (in the sense of distributions). The remaining collinear 
divergence is eliminated by subtracting distribution functions. 

3.2 Distribution function 

The matrix element for the DIS process q(p)+Y(l) — > l(P)+9(k) is given by (Q 2 = —q 2 ) 



\M\l%, 



-4C F <7 2 



2u(Q 2 - A 2 ) 



-2C F g[ 



2 u(t — A 



st 

2\2 



t 2 



\ 2 Q 2 



^Q 2 
t 2 



(3.9) 
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with s = (p + q) 2 , t = (p — k) 2 , u = (p — p) 2 . The resulting real corrections to the quark 
distribution function read (x = 1 — x) 



C F a s (Q) 



2ir 



1 + x 2 - 2£ x(l + x - £x) 9/ 

— — - 6£x 2 (2 - 3£x) 



In 



ix 2 



+ 1^? + ( " ~A X? - 2(1 - e)*^ + 3x ( *-J g)a 



f a? - 1 ' 2x 3 

2 X — 



12£x" 



2„3 



x(l — £x) 

(3.10) 



+ 6^xx(l - £x) - 6£V } @(x - £x 



The virtual gluon corrections are 



C F a s (Q) 
2vr 



5(1 -z) 2(1 -0 5 



Li 2 (e)-^ln 2 C + ln^ln(l-e)-y 



-31n£-- + 2£ln£ + 2£ . 



(3-11) 



3.3 Cutoff dependence 

Subtracting the distribution function according to eq. 



we arrive at the well-known 



IR finite result 24 



Uqq{z,£) = — < 2(1+2 ) 



ln(l 



1-z 



+ + ir^7- 42 - 6+i(1 - z) { 1 + 1 r}} <3 ' 12) 



up to terms in £ that vanish when £ — > 0. We are now interested in the size of this 
remainder. Again, it is convenient to calculate these corrections for the moments rather 
than the distribution in z. We take moments separately for the DY cross section and 
the distribution function. The upper limits in the z- and x-integrals are determined by 



(1 + Vf) 2 



1-2VI + .. 



(3.13) 



The linear dependence on A in the upper limit for the DY process reflects the smaller 
phase space in this case compared to the DIS cross section. For z — x ~ 1, a gluon with 
mass 



^(l-z) 2 <A 2 <Q 2 (l 



(3.14) 



can be emitted in the DIS process but not in DY. For the purpose of identifying leading 
power corrections, expansion of the moments for small A is sufficient. We obtain 
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W {1 \N,Z) = dzz N - 1 \w {1) ' Tea \z,0 + W il) ' viTt {z,0 



(3.15) 



C F a s {Q) 



7T 



tf(JV) + tf(JV + 2) + 2 1E - - 



In £ + Mjyy 



N 

2 +1 



iV 2 - N(V(N) + V(N + 2) + 2j E ) +3N-1 £ln£ + 0(0 



C F a s (Q) 



2vr 



^(iY) + ^(iY + 2) + 2 7E 



+ 



N - 2(V(N + 1) + V(N + 2) + 2 7E ) + 8 




N + 2 



where the £- independent terms Mdy and Mdis are not of interest presently. is the 

derivative of the logarithm of the T-function and 7^ = 0.577. . .. For the hard cross 
section the result reads 



C F a s {Q) 



7T 



M DY - Mdis + 



N 

T + 1 



£ln 2 £ + N 2 + 0(NlnN) £ln£ + C>(£) 



(3.17) 



The terms written out explicitly come entirely from the Drell-Yan process. 

Remarkably, the moments have no contribution linear on the cutoff A (no y/^-term) . 
Contrary to the what has been concluded from the double-logarithmic approximation to 
the matrix elements, the full hard cross section for the Drell-Yan process does not indicate 
'higher- twist' corrections of order Aqcd/Q from its IR cutoff dependence. According to 
Sect. 2.3 this also implies that there is no IR renormalon located at u = 1/2 and no 
ambiguity of order Aqcd/Q in the hard cross section. With the use of eq. (|2.25|) , we find 
the leading pole at u = 1 to be a double pole with residue proportional to A" from the 
£ In 2 £-term above. For A" > ln(Q/AQCD) the single pole with residue proportional to A^ 2 
is dominant. From eq. ( |2.9| ), we obtain the ambiguity (due to the Landau pole in ct s ) 



(3.18) 

The distribution function f q / q is less singular, both in £ and for large N. There is no 
double pole in the Borel transform and the expansion runs in AT£ (i.e. NAq CB /Q 2 ) 
rather than N 2 ^ as for the DY cross section. All results can be translated into z-space 
by the correspondence A" <-» (1 — z)~ l and similarly for the distribution functions in x. 

For completeness, we give the leading large- A" asymptotics of non-analytic terms in 
the expansion of W^(N, £) for small £ (This can most easily be derived from eq. ( |5.11 ) 
below.), 
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W^iN^)^ 1 UpasW) 21njVln^+ £ (iV 2 £) m ^ O • (3.19) 

7T [ ,^i m H] J 

The In £ and £ In £ terms are in agreement with the large- N limits of the corresponding 
terms in eq. ( |3.15| ). The double-logarithmic term £ In 2 £ is not seen because its coefficient 
is suppressed by l/N compared to the single logarithm. 

A similar analysis for the distribution function shows that the small-£ expansion runs 
in iV£ in this case. Thus, the leading IR corrections for large N for the hard cross section 
u> q q(N, a s (Q)) coincide with those inferred from eq. ( |3.19| ) except for subtraction of the 
IR divergent ln£. 

It is interesting that linear terms in A (a/D are absent, although the reduction of the 
phase space for real emission is linear in the cutoff. From the phase space reduction it is 
evident that any potential linear term must originate from z close to one, that is from soft 
gluons. The cancellation of such terms is possible, because the distribution W^' Tea " l (z, £), 
although finite at the endpoint of the z-integration for any finite A, is highly singular 
in the limit A — ► in the region where z is close to 1. For the same reason it is not 
legitimate to expand the distribution in A before taking the ^-integral. Expansion of the 
integrand would yield a series in £/(l — z) which is singular in the endpoint region. In 
this sense the absence of linear terms is due to a cancellation between the modification 
of phase space and the distribution over z close to the endpoint in z. 

Let us also add the following important comment: The expansion of the moments 
above assumes iV 2 £ < 1. At the same time we did not specify that N is large (compared 
to unity) and the result is applicable for small as well as large N and therefore indifferent 
to the semi-inclusive limit. Since, as a physical IR cutoff, A should scale with Aqcd, our 
conclusions are valid as long as N < Q/Aq CD . They apply to the limit N 3> 1 as long 
as one does not enter the domain of very large N, where no perturbative approximation 
is possible. This is the domain when the typical transverse momentum of the emitted 
gluon(s) is of order or below Aqcd- At this point the language of power corrections looses 
its meaning, since all such 'corrections' are of order unity and there is no leading term 
to start with. The absence of 1/Q-corrections thus holds over the entire perturbative 
domain of N. 



4 Double-logarithmic versus soft limit 

In the previous section we concluded that the hard Drell-Yan cross section in the one-loop 
approximation does not depend linearly on the cutoff and that there is no corresponding 
IR renormalon indicating a l/Q-correction. On the other hand, in Sect. 2.1 we have seen 
that such a correction arises from the same one-loop diagrams in the double-logarithmic 
approximation (DLA), that is the phase space region, where the emitted gluon is soft 
and collinear. In this section we take a closer look into this apparent difference. 

Rather than a finite gluon mass, we now use a cutoff fi on transverse momentum 
k t = uj sin 9 and energy u of the emitted gluon. We work in the soft approximation 
uj = Q(l — z)/2 <C Q. For comparison with the DLA, it is sufficient to keep only the first 
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term in square brackets of the matrix element in eq. ( |3.5| ). It can be written as 2Q 2 . 
The phase space integral eq. (|3.6|) is given by 



,2 



W^(z, fi) = S(u - M ) / S» ; Q ■ (4.1) 

Taking moments, 



l-WQ Q 2 (l~z) 2 /A 

W P)**(N tlM ) = 2 9l^M J dzz N-i J . j (1.2) 



^(1 - z) 2 - 4fc t 2 /Q2 



In the DLA one restricts fcj <C u;, since only this region gives rise to In 2 N. Then the 
/c 2 -term in the square root can be dropped and we obtain 



vr J 1 - z J fcf 

o M 2 

Adding virtual corrections 

(z N-l ^ z iV-l 

— 1) and subtracting the distribution functions, 
one obtains the standard expression, eq. (|2.5|) . From the previous section it is clear that 
a linear term in the cutoff does not arise from the DIS process and the virtual corrections 
in DY. We can therefore work with the previous expression, although the limit \i — > 
can not be taken. It is understood that the logarithmic divergences that arise in this 
limit would be cancelled by virtual corrections and subtraction of distribution functions. 
To find the dependence on the cutoff we take the derivative 

1-2/VQ 

M 2 ^ W (l),DLA (iv ) = 2 CWQ) f dz Z_^_ 

du z IT J 1 — z 



and rewrite z N ~ l = (1 — (1 — z)) N ~ l , then using the binomial theorem. This gives 

M 2^w(i),DLA (iV) ) „ 2 C F a s (Q) m (45) 

d[l A TT Q 

where the symbol £ ~' means that we keep only linear terms in [i (and ignore also ln^u 2 ). 
As anticipated one obtains linear cutoff dependence from the soft-collinear region k t <C 
u < Q. 

Let us now relax the assumption of collinearity and allow k t ~ uj. In this case (from 
eq. (H)) 

The square root can in fact not be approximated. When expanded 
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2fi- 

v/(l - z) 2 - V/Q 2 1-z ' Q 2 {T-z 



all terms are of order one in \x after integration over z. We can still use this expansion, 
provided we resum all terms of order \x. As before, we first obtain 



2 jL w mmu n = 2 gWg) yl r(i/2) 
M <fc 2 1 ' Mj * fc! r(i/2 - k) 








(- 


Q 2 ) 











(4.8) 



This expression is invalid for I = 2k. However, we do not need this case, because linear 
terms in \x originate only from 1 = 1. Thus 

,) ~ 2 M (iv-i)ig . ( , 9) 

The DLA gives precisely the first term in the sum. Recognizing the coefficients of the 
series as the Taylor-coefficients of (1 + rr) 1 / 2 , the sum of all terms gives 

^ J_ ^(i),soft (iV) ^ ^ 2 gWQ) 2 _ ^ + 1/2 = Q 

All linear dependence on /x has disappeared. Consistent with the exact calculation 
of the hard cross section, power corrections are of order N 2 fx 2 /Q 2 . As it turned out, 
corrections of order N^l/Q arise also from the regions of phase space when a soft gluon 
is emitted with large angle (k t ~ a;) and no corrections are present after adding all 
regions. Another way to state this conclusion is that corrections of order 1/Q arise 
as an artefact of the collinear approximation, which is valid to the accuracy of leading 
logarithms in N, but not to power-like accuracy. This result has important implications: 
While soft and collinear gluon emission has many universal features that allow to resum 
leading large logarithms for many processes by a single universal anomalous dimension, 
the same universality does not extend to power corrections TVAqcd/Q- Depending on 
the particular process considered the corrections of this type inferred from the soft and 
collinear region might or might not be cancelled by those from other regions. 



5 Soft factorization in the single-chain approxima- 
tion 

The DY cross section calculated to first order in a s with an explicit IR regulator does not 
contain IR contributions linear in the cutoff. In this section we translate this result into 
the language of large-order behaviour in perturbation theory. This serves as a check that 
the absence of contributions of order 1/Q is consistent with perturbative factorization 
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0, [3], |I4" |, and elucidates which of the options A, B, C, of Sect. 2.2 is realized. We have 
seen that for the cancellation of contributions proportional to N/Q it is necessary to 
account exactly for soft gluon emission at all angles. The large logarithms in TV beyond 
the DLA also exponentiate || |3], |14| and in fact correct for soft gluon emission at large 



angles, but it is only after resummation of all subdominant logarithms that the result 
is adequate to exhibit the cancellation of all IR contributions of order 1/Q. Hence, to 
clarify the situation the anomalous dimension functions in the exponent of the resummed 
cross section have to be calculated to all orders in perturbation theory, in agreement with 
general considerations in Sect. 2.2. 

Since we do not expect the DIS subprocess to cause power corrections of order 1/Q, 
we restrict ourselves to the contributions from the DY process, contained in the functions 



A and C in eq. (|2.2j ). Being interested only in the soft gluon region, we recall that the 
emission of soft gluons from the incoming quark and anti-quark can be accounted for by 
eikonal phases equivalent to Wilson lines along the classical trajectories of the partons 



14| . This is of particular interest in the present context, since the Wilson line operators 
take into account all soft contributions, including the non-collinear ones. In view of the 
previous discussion, we would therefore not expect any 1 / Q-corrections to emerge in this 
approach. 

5.1 Wilson lines 



We follow the treatment of Korchemsky and M archesini ]T4 ] , generalizing it to our par- 



ticular class of higher order contributions. For more details, we refer the reader to the 



second reference of [14 . 



The matrix element for emission of soft partons is given by (?2|T£/dy(0)|0), where 

Udy(x) = -Pexp ^ig s J ds p^A /1 (p 2 s + x)^j Pexp ^-ig s J dsp^A^pis + x) j , (5.1) 

is the product of Wilson line operators describing the annihilation of an on-shell quark 
and anti-quark with momenta p\ and p 2 at the space-time point x. Up to corrections 
that vanish as z — > 1, the partonic Drell-Yan cross section is given by 

W(z,a s (Q)) = H DY (a s (Q))W DY (z,a s (Q)) . (5.2) 

Hjyy = 1 + 0(a s (Q)) is a short- distance dominated function, independent of z. Wby is 
the square of the matrix element, summed over all final states: 



CO 



W DY (z,a s (Q)) = ^ J J£ e <»>Q(i-*)/2 (0|T C/J Y (y) T [/ DY (0)|0) (5.3) 

— oo 

The Fourier transform is taken with respect to the energy of soft partons and y = (yo, 0). 

The crucial observation is that the Wilson line depends only on the ratio (/j,N) / (QNq) 
(taking moments of Wr>y(z, a s (Q)), where \x is a cutoff separating soft and hard emission 
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(the renormalization scale for the Wilson line) and N Q is a suitable constant. Hence the 
iV-dependence of the Wilson lines can be obtained from their //-dependence, which is 
given by the renormalization group equation 



2 d a ,_ n d \ , Tir ( u 2 N 2 _ , A \ ^ , u 2 N 2 



»W 2+(3{as) d^s) ^^W^r^) = W«-)ln^ + r DY (a s ). (5.4) 



Here r cusp (a s ) is the universal cusp anomalous dimension of the Wilson line |[26|| . The 
general solution to eq. ( |5.4| ) is given by 



I i, 2 N 2 \ r dk 2 

lnW DY (^,a s (u)\=X(a s (QN /N))+j 2* 



Q 2 V 2 /7V 2 * 



x cusp 

(a s (k t ) In -2—2 + T DY (a s (k t )) 

(5.5) 



The integral is a particular solution of the inhomogeneous eq. (|5.4j) and X(a s (QNo/N)) 
denotes the general solution of the corresponding homogeneous equation, an arbitrary 
function of the running coupling a s (QN /N). Putting \i = QN /N, one identifies 
X(a s (QN /N)) = In Wdy(1> a s (QN /N)). The inhomogeneous term can be rewritten 
(identically) in a more familiar form and we obtain for the partonic Drell-Yan cross 
section 



l-No/N 

dz 



W(N,a s (Q)) = H DY (a s (Q))W DY (l,a s (QN /N)) exp I 2 j - 



Q r dk 2 \ 

J -^r cusp (a s (k t )) + T DY (a s ((l-z)Q)\ + OQn N/N) . (5.6) 

(l-z) 2 Q 2 * ' 

We have restored the hard part H-D Y (a s (Q)) and set fi — Q. This equation differs in 
form from the standard resummed cross section eq. ( |2.1| ) with the exponent eq. ( |2.2| ) 
written with the replacement eq. (|2.8|) only by the presence of the initial condition 
Wdy(1, a s (QN /N)). Its expansion in a s produces subdominant logarithms a h \n h ~ 1 N. 
To fully conform with the conventional expression for the resummed cross section, the 
initial condition can be absorbed into the following redefinitons of Hd Y and T-£, Y : 



H DY (a s (Q)) — » H DY (a s (Q))W DY (l,a s (Q))) 

r D y(a s ) — ► T DY (a s ) - (3(a s )J- lnW DY {l,a s ) (5.7) 

da s 

The redefined Tdy starts at order a 2 . It does not affect resummation of large logarithms 
in N to next-to-leading accuracy a s (Q) k ln fc N. The argument of the coupling in Tdy is 
naturally (l — z)Q and not \/l — zQ as in the function B in eq. ( |2.2|) . It is for this reason 
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that we have reintroduced the function C there. If we absorbed C into a redefintion of 



A and B as in [17]], A would no longer coincide with the universal function r cusp starting 
from order a^. 

In the remainder of this section we reproduce the soft factorization formula (|5.2j ) to 
all orders in the strong coupling, but in the single-chain approximation. We do so first by 
taking the large- N limit of the result obtained in Sect. 3, and then by direct calculation 
of the Wilson line. We will see that 



- IR renormalons in the Wilson line are in exact correspondence with the IR renor- 
malons in the full cross section in the large- N limit. This excludes the scenario C in 
Sect. 2.2 which would have involved the remainder R, which is suppressed by 1/N to 
any finite order in perturbation theory and not captured by eq. ( |5.2| ). 

- the Wilson line satisfies the correct RG equation, with the anomalous dimensions 
r c us P (a,s(Q)), T^y(o: s (Q)) being entire functions in the MS scheme. IR renormalons enter 
exclusively via the initial condition for the evolution of the Wilson line, which by com- 
parison with the calculation in the single-chain approximation requires non-perturbative 
corrections only at the level of 1/Q 2 - This verifies the absence of l/Q-corrections to the 
DY cross section within the soft factorization technique. 

- the standard exponentiated cross section chooses a particular solution to the evo- 
lution equation that requires a redefinition of both the anomalous dimension functions 
and the boundary condition in such a way that both of them contain an IR renormalon 
(pretending 1 / Q-ambiguities) that was not present in the initial formulation and which 
cancels in the sum. This cancellation selects scenario B from Sect. 2.2. 



5.2 Large- N limit of the DY cross section 

In this subsection we calculate the cross section W(z, a s (Q)) on the l.h.s. of eq. ( ^2|) 
as the large- N limit of the result obtained in Sect. 3. Instead of calculating the series of 
higher orders in a s , it is more convenient to study the Borel transform, eq. ( p. 11 ), which 
can be obtained by Mellin transformation, eq. ( |2.25| ), from the one-loop result with finite 
gluon mass, as given in Sect. 3.1. The virtual corrections are iV-independent and of no 
importance for what follows. For the contribution from real emission we obtain 



B 



W 



real 



(N,u) 



Cf sin(7ru) 



3 5u/3 



7T 



ITU 



x z 



N-l 



In 



v 



z ■ 



£z-y/{z + Szy-4£z 



(5i 



neglecting all terms that do not contribute to the leading asymptotics at iV — >• oo. After 
repeated substitution of variables, the integral is 
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B 



W 



real 



(N,u) 



C F sin(Tra) ^ u/3 

TT 7TU 



1 + 



■ '2/V / In ; ' v ~ (l-z)- 1 -" 



1- 

u 

5(1 - 2u, 2JV + 2u + 1)2^1 Q-u,l-«;iV+|;x 
+ (1 - x)£(2 - 2u, 2iV + 2u) 2-F1 (J - «, 1 - «; iV + |; x 



(5.9) 



where B(x,y) = T(x)T(y) /T(x + y) and 2-^1(0, c; x) denotes the hypergeometric func- 
tion. Because of the arguments of B, the second term in square brackets contributes 
only at relative order 1/N can be dropped. The first term can be integrated using 



The result is 



N + 1/2 d 
u(N + u) dx 



(1-x) 



—u 77 
2-fl 



u, —u; N + -; x 



= (l-x)- l - u 2 F 1 [^-u,l-u;N + ^;x) . 



(5.10) 



B 



W 



real 



(N,u) 



Cp smjiru) 5u/3 r(-u) 2 r(N + uf 



TT TTU 



T(N)< 



N^l Cf_ 5u /3 J_ T ( l ~ u ) 



TT U 2 T(1 + U) 



2u 



(5.11) 



Close to the poles at integer u = m the result is accurate to the dominant power of N, 
i.e. N 2m , in the residue, neglecting corrections of order (In N)/N. The double pole at 
u = is cancelled by the virtual corrections. Poles at half-intergers do not arise. 

A similar analysis for the DIS distribution function shows that the iV-dependence of 
the Borel transform in the large- N limit is given by N u in accordance with the observation 
in Sect. 3.3 that the small-£ expansion runs in iV£ in this case. 



5.3 Calculation of the Wilson line 

Next, we obtain the Wilson line W^dy- The Wilson line operators need renormalization 
and we choose MS subtractions. This implies in particular that all integrals without 
a scale are set to zero. Since p\ = pi = 0, all virtual corrections to Wby vanish and 
non-zero contributions arise only from contraction of gluon fields with different time 
ordering. The leading-order expression is then given by 



w«/ \ n 2<3 (f?e»*X f d d k 2 f Q ,\ 2 Pl ■ p 2 

W^(z,a s ) = C F g s — I I J —— d 2TT5{k)5[k --{l-z 




pi ■ kp 2 ■ k 



t-WWD <5 - 12) 
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Figure 2: (a) Cut diagrams with fermion loops, (b) No n- vanishing contribution with 
no cut fermion loop. The crosses denote counterterm insertions. 



Note that we have inserted the factor e 7E /(47r) in front of fi 2 . With this choice of 
renormalization scale /i, MS corresponds to subtraction of poles in e (d = 4 — 2e). 

To compute diagrams with fermion loop insertions into the gluon line, we follow 
Appendix A of [25|. Each fermion loop in a diagram like in Fig. is accompanied by a 
factor 



«s(m) 



0& 



L{e) 



Lie) 



6r(l + e)r(2-e) 2 
r(4 - 2e) 



(5.13) 



and /3g = Nf/(6ir). As in lowest order, the only non- vanishing type of diagram is that 
shown in Fig. ^a and the symmetric one. The cuts in this diagram proportional to 5(k 2 ) 
all vanish except for the diagram with no fermion loop, since with k 2 = the fermion 
loop is given by an integral without scale. The diagram with no loop will be added 
later. The non-zero cuts are those proportional to 0(/c 2 ), labelled (b), when a fermion 
loop is cut. The corresponding imaginary part of the effective gluon propagator with an 
arbitrary number of fermion loops is0 



oo 
n=l 




L(e) n 



smir(-ne) 



(5.14) 



The k^k u -terms can be dropped due to gauge invariance (they are cancelled by diagrams 
not shown). A short calculation gives for the sum of all diagrams with cut fermion loops 

8 A pedestrian way of deriving this is to sum and square the diagrams with one cut bubble literally. 
With the correct phases, the sum of all interference terms in order a™ combine as 



n-l 



„/>6(l-e)r(2-e) v-^ 



T(4 - 2e) 



k=0 



l))ne = —L(e) sin(— nire) . 
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C F a s 



7T 1 — Z 



E< 

n=l 



(n + l)(-e) 



n+l 



G(-e,-(n + l)e), 



(5.15) 



-(n+l)e 



T(l/2 - (n + l)e)r(l + ne) 



When we add the counterterms for the fermion loops, a non-vanishing contribution 
requires again that at least one fermion loop is cut, except when all fermion loops are 
replaced by counterterms as in Fig. [|b. The second contribution will be added shortly. 



Accounting for the first one, the previous expression is modified to [25 ] 



C F a s 



1 



(-l) fc 



■n 



7T 



oo n— 1 



G(-e,-(n + l-k)e) 



(5.16) 



Next we note that the remaining contribution from the lowest order diagram together 
with the ones where all fermion loops are replaced by counterterms is simply 



1 - a s /?5/e 



(5.17) 



with W§ Y (z, a s ) as in eq. (|5.12j ). Adding this term amounts to extending the sum over 
n to and the sum over k to n in eq. ( |5.16| ).P| 

Eq. ( |5.16D contains poles a™/e n+1 , but an additional 1/e is in fact present as seen by 
introducing -(—distributions. At this point it is convenient to pass to moments, before 
we take the final overall subtractions. The integral over z is trivial and we obtain 



W DY (N, a s 



oo n 

2^ e «?(/tfr E 



7T 



n=0 



o (-e) n+2 (n + l - W 



H(-e,-{n + l 



k)e), 
(5.18) 



9 Let us add an inte resti ng side-remark: For IR safe quantities like the hard cross section in DY, the 
contribution from eq. ( 5.17 ) is necessary to cancel all IR poles in e. However, in this case it is possible 



to take the Borel tr ansfor m of the series, set e to zero at this stage of the calculation and let u act as 
regulator. Then eq. (5.17) is proportional to exp(— u/e), which can be set to zero as e — > 0, provided we 



assume u to be (small and) positive. Only the Borel transform of a term like eq. ( p,16| ) is left, which 
apparently does not contain the lowes t ord er contribution in a s . However, due to an additional l/u, 
the Borel transform of a term like eq. (5.1£) is in fact of order one close to u = 0, so that re-expansion 
of the Borel integral in a s produces an order a s correction, with coefficient equal to the lowest order 
correction. It is not difficult to derive eq. ( 2.25 ) in this way. Alternatively, one can follow this section, 
put e = from the very beginning and use 



(2tt)4 



= / d\ 2 



d 3 k 



to obtain the Borel transform in the Mellin representation of eq. ( [2.25 ) 
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ff(e ' s) = [L( - 6)] [W^) 2T(N + 2s)T(l-s + e) - (5 ' 19) 
Note the signs of the arguments of the function H. We introduce expansions 

OO CO 

H(e,s) =Y,HM j , H s {e) = £ hf* . (5.20) 

j'=0 fc=0 

The functions ifj are finite for e — > by construction. For later use we collect the 
expressions 

Q 2 -5/3V r(JV)r(i + fl ) 



H(0,s) = i^e 



/i 2 y 2T(N + 2a)r(l - s) 

■a- 1 \ _ i r(4 + 26) 

^ o(ej " 2 6r(i- e )[r(2 + e )pr(i + e) (5 ' 2ij 

#l(e) = H (e)|ilnL(-e) + ln^^-2*(JV) + *(l + e)| . 



To subtract the overall divergences from eq. ( |5.18| ), H(—e, —(n + l — k)e) is first expanded 
in its second variable 



7r 

n=0 



2 ^^ a „ ( ^ r j^ wn+2) _^ (1)) _^ F0( _ e) (5 22) 



To arrive at the second line we use 



fc=0 



(-l)'7(n + l)(*(n + 2)-*(l)) j = 
(_!)»/(„ + 1) j = i 

j = 2,...,n + l 

n! j = n + 2 

(5.23) 

Minimal subtraction of the pole part in e gives the final expression for the MS- 
renormalized Wilson line with an infinite number of fermion loops: 

W DY (N,a s ) = 2^f;<(/3 / r{t^(*(n + 2)-*(l))/ l M 2 

+ ^T /l I 1 +i + n!F «+2(°)} ( 5 - 24 ) 
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Note that the result contains two terms related to the 'anomalous dimension functions' 
H and Hi, that appeared as coefficients of the l/e n+2 and l/e n+1 pieces before, and a 
'finite term' with coefficient n\. It is instructive to take the Borel transform of the above 
series. We then obtain 



B[W DY ] (N,u) 



7T 



H(0, -u) - flo(0) + uHi{u) 1 d f ^xH {u) - H (xu) 



u du 



dx- 



x(l 



■ x) 



(5.25) 



Y<n=o h]i'u n /n\ is essentially the Borel transform of Hj(e), cf. eq. ( |5.20|) . 

from eq. (|5.21 



-u 



we find that this term above 



where Hj(u) 

With the explicit expression for H(0, 
coincides with the large- iV limit of the Borel transform of the full Drell-Yan cross section 
in eq. ( |5.11|) , up to corrections 1/N that have been neglected in both cases. There it 
was understood that the double pole at u = would be cancelled by virtual correc- 
tions and distribution functions, but eq. ( |5.11|) provided the leading power corrections 
(iVA Q cD/Q) 2fc (k = 1,...). In the present case the divergences are subtracted by the 
functions H and Hi. These additional terms are analytic in u (as expected for anomalous 
dimensions in minimal subtraction schemes), so that indeed the Wilson line operators 
incorporate the leading power corrections (with the largest positive power of N) to the 
hard Drell-Yan cross section. In particular, no power correction of order 1/Q appears. 
Another important observation is that only the function Hq is related to the most singu- 
lar pole in e, which comes from soft and collinear emission. While H Q will therefore be 
related to the universal cusp (eikonal) anomalous dimension, all conclusions on power 
corrections follow from H (0, —u) which is process-specific (that is, depends on the Wilson 
contour and not its cusps alone). 



5.4 Evolution equation and resummation 



Since eq. ( |5.24|) is exact in the single-chain approximation, it must satisfy the RG equa- 



tion Q5.4D to the same accuracy, that is with /3(a s ) = Poa^ and ln[l + Wby] equivalent 
to Wjjyi see the discussion in Sect. 2.3. (The unity under the logarithm is the tree-level 
contribution W^ Y — 1-) Using the explicit expression in eq. ( |5.25| ) we get 




1 r°° 
— / due u/ ^ oas) 



0o Jo 



(l 



d 



£arfli[*(" + 2)-*(l)] 

n=0 

J_ f 1 dx H (a s ) - H (xa s 
a s Jo I- x 



B[W DY ](N,u) 



(5.26) 



where we have restored the full non-abelian fl = —(11 — 2/3A^)/(4-7r) and introduced 

a s = -0O&S ■ ( 5 - 27 ) 
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In the large- iV limit we expect to find a term logarithmic in N, multiplied by the cusp 
anomalous dimension. From the explicit expression for Hi we get indeed 



2C F a s Cfol s 
-Hi(a s ) ■■ 



T(4 + 2a s 



7T 



vr 6r(l-a s )[r(2 + a s )] 2 r(l + < 
with N = e~ lE and omitting terms without IniV. Thus 

C F a s r(4 + 2a. 



• 2 In 



(5.28) 



cusp 



C F a s 

2 H (a s 

vr 

C F a s 



vr 6r(l -a s )[r(2 + a s )] 2 r(l + a s 



VT 



1 - -a. 



-a 2 . + 2.0708ai? - 1.0932a 4 e 



(5.29) 



in agreement with the known two-loop result (keeping the leading Nf-term only) [IR |26 
We verify coincidence with the direct calculation in Appendix A. 

The remaining terms in eq. (|5.26| ) (up to 1/iV-corrections, which are consistently 
neglected) give Tdy^s): 



2 C F a s ^ 


H (a s 
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vr 


) 




a s 


C F a s 




vr 




a, 


vr 1 




24 j 



lnL(-a s ) - ^(1 + a. 

0.08547o? - 0.1537a? + 0.4951at + . . . 



n=0 



(5.30) 



Note that r DY (a< j ) starts from second order in a s fl4}| . Eq. (|5.25|) can now be written as 



B[W DY ](N,u) = B[W rcal ]{N,u) 



0_ 

Q 2 



1 



u- 



-B[r cusp ](u) 



1 + 2nln 



1 



5[r DY ](«) 



(5.31) 



The first term on the right hand side gives the unrenormalized Borel transform, which co- 
incides with I?[W real ], eq. (|5.11 ), and the two other terms are the subtractions, necessary 
for finiteness at u = 0. 

The anomalous dimensions, eq. ( |5.29 ) and ( p. 30 ), have finite radius of convergence 
(in the MS scheme) and their Borel transforms are entire functions. At u = 1/2 they 
are non-zero. We conclude that none of the terms in the evolution equation for the 



Wilson line, eq. ( 5.4|) , contains renormalons, and neither does the particular solution 
corresponding to the integral in eq. ( |5.5|) . (Let us insist again that one must stay in the 
perturbative domain iV < Q/Aqcd)- In the present formalism, all IR renormalons to 
resummed cross sections are introduced through the initial condition (or solution to the 
homogeneous equation) for the evolution of the Wilson line. The initial condition can be 
determined by comparison with exact calculation (in the single-chain approximation), 



1 



N 



B[X](N,u) = B[W™ l ](N,u) - -j ^— ) [B[T cusp ](u) + uB[T DY ](u)] , 



2u 



(5.32) 



2s 



and consequently, the renormalons in the initial condition are exactly those in eq. ( |5.25[ ). 
We find once more that the dominant IR renormalon requires non-perturbative con- 
tributiuons only at order N 2 Aq CD /Q 2 . All terms in eq. ( |5.32| ) have the same N- 
dependence proportional to N 2u . This verifies that the initial condition is a function 
only of a s (QN /N). Indeed, using eq. ( |2.12| ) we get 



B[X(a s {QN /N)](u) 



N J 



2 a 



B[X(a s (Q))](u) 



(5.33) 



It remains to clarify the origin of apparent linear power corrections to the resummed 
cross section in the form of eq. ( |2.2|) . The difference with eq. ([5.5|) lies in the choice of 
the inhomogeneous term in the solution to the evolution equation as 



f V 2 N 2 N 

lnT Mw as(/i) 



X'(a s (QN /N)) -2 dz 



} , z N ~ l - 1 


r f 


2 / dz 
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J 1-z 







L (i-2) 2 <3 2 



+ r DY (a s ((i-^)Q) 



dkf „ 
k 2 usp 



{a s {h) 
(5.34) 



A simple calculation (as in Sect. 2.2) yields for the difference of Borel transforms of the 
integrals in eq. ( |5.5|) and eq. ( |5.34| ) 



N \ 2u / \ 1 

—J (iV - 2 «r(l - 2u) - l) ^2 [^[ r cu SP ] W + uBFdyKu)] . (5.35) 



This difference has to be compensated by a change in the solution of the homogeneous 
RG equation, which now becomes, 



B[X'](N,u) = B[W rcal ](N,u) 



N 



2h 



N ~ 2u r(i 



2u) [B [r cusp ] (u) + uB [Tjyy] (u)] . 

(5.36) 



Comparing with eq. ( |5.32j ), we see that this change involves singular terms at u = 1/2. 
Because of the r(l — 2u), both, the Borel transform of X' and of the integral in eq. (|5.34|) 
have IR renormalon poles at u = 1/2 (corresponding to 1 /Q) with residue proportional to 
N, but opposite sign, so that in the sum the singularity cancels. In the Borel transform of 
the integral in eq. ( |5.34| ) the pole is present through the functions Ai and A 3 , introduced 
in Sect. 2.2, and is cancelled in this representation by an explicit pole in X' . If, to return 
to the standard representation, the initial condition contribution X' is eliminated by 
the redefinitions of eq. ( |5.7| ), this cancellation persists in the exponent. Indeed, for the 
functions A and C in eq. (|2.2j) we get 



B[A](u) 



B[r cusp ](u) 
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-B[C](u) = B[T DY }(u)+u^ T{1 l 2u f [X}(N,u) (5.37) 



AT \ 2u /V"2it | 

„ c ^.. ]Wa) (_JL) „,„. 



so that 

2B[A](u) - uB[C](u) = r(i ° 2u) u 2 (^J B[^ cal ](iV, u ). (5.38) 

Due to the Gamma-function in the denominator, this combination has a zero at u = 1/2, 
independent of the particular form of B[W TeaX ](N, u) (unless it is singular, which is 
not the case, at least in the single-chain approximation, see eq. (|5.25|) ). Note that all 
contributions from the anomalous dimensions r cusp and T^yy disappear. 

We conclude that the cancellation of scenario B in Sect. 2.2 does occur and no indi- 
cation for physical 'higher twist' corrections of order 1/Q remains. The cancellation of 
such corrections with the remainder R as in mechanism C could only have taken place, 
if the soft approximation had accounted correctly for the leading- N behaviour only to 
logarithmic, but not to power-like accuracy. This is not the case. As noted before (and 
expected), W^y computed in the eikonal approximation (Wilson line) reproduces those 
power corrections to the Drell-Yan hard cross section with the largest power in N. This 
suggests that the nonperturbative corrections indicated by renormalons in the Wilson 
line should exponentiate |5|]. 

6 Discussion and Summary 

Let us summarize our investigation of power corrections (alias infrared renormalons) in 
Drell-Yan production. Our starting observation is that the moments of the resummed 
hard cross section in the large- N limit can be written in two different ways as 

2 (l-z) 2 

dk 2 




lnu qq -(N,a s (Q)) = - / dz - — \2 / -± T cnsp (a s {k t )) + B(a s {VT^Q)) 



+ C(a e ((l-z)Q)\+0(l) (6.1) 



or 



l-e-TB/JV Q 2 {l-z) 



r if/" dk 2 

lnu qq {N,a s {Q)) = J dz Y~\ 2 J ~W V ™v^s(k t )) + B{a 3 (VT^Q)) 



Q 2 (l-z) 

+ C(a s ((l-z)Q)\+0(l). (6.2) 
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Both forms are equivalent for the purpose of summing large logarithms in N, but have 
apparently different properties as far as sensitivity to the IR Landau pole of the coupling 
is concerned. When the 'anomalous dimensions' r cusp , B etc. are truncated at finite 
order, the second form has no ambiguity unless iV > Q/Aqcd, whereas such an ambiguity 
is present in the first version for any N. This apparent conflict is resolved by the analysis 
of the exponent in Sect. 2.2, which concludes that the correct IR renormalon structure 
of the exponent can be obtained only, if the anomalous dimensions are also considered 
to all orders. With respect to eqs. ( |6.1| ), ( |6.2|) this means that IR renormalons that are 
manifest in terms of the integrals over z and k t (the functions A, in Sect. 2.2) in one 
form can be hidden in the large-order behaviour of the anomalous dimension functions in 
the other. It is with this observation that we start to differ from previous works, which 
implicitly assumed that low-order approximations to the anomalous dimensions would 
be sufficient Q or the functions B and C, which account for gluon emission that is not 
both, soft and collinear, could be neglected || |T0|j . 

Already from this conclusion it is evident that, even in the large- N limit the search 
for power corrections is a problem quite different from resummation of In N, which can 
be completed up to a desired order with finite-order approximations to the anomalous 
dimensions. Although the techniques of soft factorization developed for the purpose of 
resummation are still effective in isolating those power corrections (NAqcd/Q) 21 * with 
the largest power of N, this task requires application of these techniques to all orders in 
perturbation theory. 

We have argued that taking into account these considerations, all corrections of order 
1/Q cancel in the Drell-Yan cross section and the remaining IR renormalons indicate 
parametrically smaller power corrections of order N 2 Aq CD /Q 2 . The evidence for this 
cancellation has been collected from two sources: A calculation of the Drell-Yan cross 
section for any N with finite gluon mass as IR regulator (which by eq. ( |2.25| ) can be 
considered also as a certain approximation to large orders in perturbation theory) does 
not reveal any IR sensitivity proportional to X/Q, even though the phase space in the 
presence of finite gluon mass is reduced linearly in X/Q. The same result has been 
obtained from implementing the factorization of soft gluon emission in terms of Wilson 
lines |14]] in the approximation of a single chain of vacuum polarization insertions. This 
approximation may be considered as a model (exact in the unphysical limit of large 
number of fermions) for the large-order behaviour of the anomalous dimension functions. 
In this language the absence of l/Q-power corrections seems to be most transparent. 
Renormalizing the Wilson line by minimal subtraction, we find that the anomalous 
dimensions that enter the renormalization group equation for the Wilson line are analytic 
for small a s and do not contain IR renormalons. These enter the solution of the evolution 
equation exclusively through the initial (boundary) condition, which, in general, contains 
both, perturbative and non-perturbative contributions. As in the exact calculation of 
the hard Drell-Yan cross section, the power corrections deduced from IR renormalons are 
of order N 2 Aq CD /Q 2 . It is only when this initial condition is absorbed into a redefinition 
of the functions B and C that IR renormalons enter the exponent. These redefintions are 
such that the power corrections of order 1/Q found from eq. ( |S.lp [||, Q are eliminated 
by an exact zero in the Borel transforms of the corresponding anomalous dimensions. 
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In Sect. 4 we have seen that the identification of power corrections through IR renor- 
malons or cutoff dependence relies on an exact treatment of soft gluon emission also at 
large angles, while the collinear approximation can introduce spurious cutoff dependence. 
It is for this reason that the leading (in N) power corrections can not be described by 
the universal cusp anomalous dimension. In terms of Wilson lines this is seen explic- 
itly by the fact that IR renormalons enter through the initial condition and not as the 
coefficients of the leading soft and collinear pole in e. We also note that the result on 
power corrections is not specific to the semi-inclusive limit and applies to all N in the 
perturbative domain N < Q/Aq CD . The analysis of Wilson lines suggests however that, 
similarly to logarithms in N, the power corrections leading in N can also exponentiate. 

All explicit calculations of large-order behaviour have been done in the approximation 
of a single 'renormalon chain'. If we consider this approximation as an approximation 
for In uj q q rather than the hard cross scetion u qq - as suggested by exponentiation of large 
logarithms beyond this approximation, then as usual the contributions from several 
chains are bound to modify the precise numerical coefficient of the iV 2 AQ CD /Q 2 -terms. 
To confirm or disprove the absence of 1/Q contributions to the Drell-Yan process at 
the level of two gluon emission would again require an exact treatment of soft gluon 
emission that goes beyond the usual approximations in treating parton cascades. From 
the fact that even in the single-chain approximation the IR renormalon structure is 
sensitive to the large-order behaviour of the anomalous dimensions, we expect this to 
be a quite challanging problem. Without any further physical evidence to the opposite 
we find it reasonable to assume that as in other cases (with operator product expansion 
as a check) || multiple chains and more complex diagrams will not modify the power 
behaviour itself. Still, we recall that the Drell-Yan process may receive power corrections, 
which can not be seen through IR renormalons, and whose size is not accessible with 
present theoretical tools. 

As pointed out in [Q, the question of linear power corrections is of particular 
importance due to their potential universality, which would allow to treat such power 
corrections to processes without operator product expansion on a phenomenological level 
akin to matrix elements of higher dimensional operators in cases with operator product 
expansion.^ From this point of view it is interesting to compare the result on power 
corrections for the Drell-Yan process with those for hadronic event shape variables in 
e + e _ -annihilation |2T], % H 0- In particular, from the point of view of summing large 
logarithms in N to next-to-leading accuracy, the thrust distribution is closely related to 
the distribution in z in the Drell-Yan process, ou qq -(z, a s (Q)), since the large logarithms 
can be accounted for by the same jet mass distribution |27|, ES|. (For a discussion in the 
present context, see |], |10|.) On the other hand, if we compute the thrust distribution 
with finite gluon mass to order a s , we obtain for the averages 

10 The statement of universality implies that the universality of IR renormalons is also preserved 
non-perturbatively. This is not always the case. For example, the 'binding energy' A that appears in 
heavy quark expansions has a universal IR renormalon ambiguity, but its numerical value (after proper 
definition) depends on the particular hadron under consideration. Similarly, one would expect the actual 
value of corrections of order 1/Q 2 in the Drell-Yan process to depend on the particular hadron target. 
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(1-T) ~ -8G^A 

7T Q 

(T^ 1 ) ~ 8G(iV-l)^^A (6.3) 

7T Q 

where G = Efc=o(-l)V(2fc + l) 2 = 0.91596. . . is Catalan's constant.^ Contrary to the 
Drell-Yan process, large 1 / Q-corrrections are present for event shapes, in particular for 
the moments of the thrust-distribution. Some caution is necessary in interpreting the 
results obtained with finite gluon mass in these cases, since the correspondence with 
IR renormalons, provided by eq. ( |2.25| ) for the DY process, does not hold, because the 
phase space of the quark pair from gluon splitting in the single chain approximation 



is not integrated unweighted. We refer to the very recent publication p9| , where the 
relevance of four-parton contributions to l/Q-terms is discussed. The difference in 1/Q- 
corrections to the DY cross section and thrust reinforces that these corrections can only 
be identified by a treatment of soft gluon emission at all angles to all orders. 

One may speculate on a deeper physical reason for this difference. One possibility is 
that IR safety of event shape variables is in a certain sense enforced 'by hand'. Depending 
on the construction, this can enhance or suppress IR regions in an almost arbitrary way 
PTJ. On the other hand, for the DY process one might be more confident that IR 
renormalons probe some aspect of the IR dynamics of QCD that eventually might find 
a rigorous treatment just as higher twist matrix elements in deep inelastic scattering. In 
this respect it is somewhat comforting that corrections to the DY process appear to be 
of order 1/Q 2 . 

The theoretical discussion has been performed in terms of moments. To obtain a 
prediction for the resummed DY cross section, the inverse Mellin transform has to be 
taken. Let us note first that there is in fact no compelling reason to evaluate the integrals 
with running coupling in eqs. ( |6.1| ) and ( |6.2|) exactly as long as, as always in practice, 
only finite-order approximations to the anomalous dimensions are available. To any 
definite accuracy in the summation of logarithms of N, one may expand a s (k t ) etc. to 
the required order and perform the integrals term by term. Then the Landau pole in the 
perturbative running coupling never poses a difficulty. This reflects our previous assertion 
that the question of power corrections necessarily deals with the all-order behaviour of 
the anomalous dimensions. 



11 The discrepancy between the coefficient — 8G and the value —4 reported in arises entirely from 
the normalization of thrust. In the region where x\ + X2 xi being the energy fractions of quark and 
antiquark) is close to its maximum value, we have, for finite £ = A 2 /Q 2 , 

T = max - 2 max ( Xl ' X2 ) 



E< Xi + X 2 + y/{2-X! -x 2 ) 2 -4£ 



One can not set £ = in the square root, since in the region of interest (2 — X2 — X2) 2 is of order £. One 
can see that the difference is a -^-effect, since with the above normalization 1 — T = in the two-jet 
region x\ + X2 — 2(1 — while with £ = in the expression for T, one would obtain 1 — T = a/J. 
See also Appendix B. 
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Nevertheless the exact integrals over the running coupling are often kept and the de- 
pendence on regulating the Landau pole is then interpreted as the ultimate accuracy that 
can be achieved without deeper understanding of non-perturbative corrections. From our 
analysis we conclude that such a treatment has to be interpreted with care, since the 
failure to include the anomalous dimensions to all orders may pretend artificially large 
regularization dependence. From this point of view, the construction of the inverse Mellin 
transform starting from eq. ( |6.2| ) is preferred compared to eq. ( |6.1[ ), which was used in 
[|J, since spurious large corrections of order NAqcd/Q do not appear. We find confir- 
mation for this suggestion from the analysis of |5t| , where comparison of the resummed 
cross section based on eq. ( |6.1|) [§], with data requires to fit rather large corrections 
proportional to NAqqd/Q, whose main effect is to reduce the effect of resummation. 

The phenomenological implications could pertain to processes other than Drell-Yan. 
For instance, a resummation similar to the Drell-Yan case has also been completed for 
the top production cross section (where in addition large corrections from gluon-gluon 
fusion have to be considered) |HJ]. It might be worthwhile to investigate the rather large 
cutoff dependence reported in the light of the results presented in this paper, especially 
since the cutoff was chosen much larger than Aqcd- An improvement of the theoretical 
Standard Model prediction is crucial to unravel potential anomalous couplings of the top 
quark in its production cross section. 

We conclude that further explorations of power corrections, both theoretical and 
phenomenological, are expected to provide further insight into the workings of QCD 
and will help understanding the limitations of perturbative QCD in describing data for 
hadronic observables. 
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Appendix A 

In this appendix we sketch a simple derivation of the cusp anomalous dimension to all 
orders in the l/A/-approximation (single-chain approximation). In lowest order, the 
cusp anomalous dimension is given by the 1/e-pole in the diagram shown in Fig. |3|. In 
higher orders we adopt dimensional regularization with minimal subtraction. As shown 
]l9f the series in the 1/A^-approximation can be extracted from the large- A behaviour 



m 



of the diagram, when evaluated with a gluon of mass A. A straightforward calculation 
gives 



V 2, / d d k pi-pa C F a s r(l + e) f X 2 \ 
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Figure 3: Diagram for the cusp anomalous dimension. The double lines denote a Wilson 
line with cusp (alternatively, they can be interpreted as eikonal propagators). 



where cosh 7 = (pi ■ V<i) I vV\!p\ is the Minkowskian cusp angle [p6| . The cusp anomalous 
dimension can be read off the coefficient of 1/e (A 2 //i 2 ) -€ ||19[ (Appendix A) and we 
obtain (adding self-energy contributions, which result in 700th 7 — > 700th 7 — 1) 



r, 



cusp 



r(4 - 2(3 a s 



vr 6r(2 - P a s ) 2 T(l - p a s )T(l + f3 a s 



(700th 7 — 1) 



C F 



as 



7T 



108 



C F N 



+ ... 



(7coth7- 1),(A.2) 



where (3q = Nf / (6tt) . The second order coefficient is in agreement with the Aj-term of 



the complete two-loop result [15, 2c]. The above series gives the term with the largest 
power of Nf to all orders. 



Appendix B 

The different power corrections to the Drell-Yan cross section and averages of thrust can 
be reduced to a few elementary integrals. First we write for the moments 

(X N - 1 ) = (1)-(N-1)(1-X) + ... , (B.l) 

where X = T (thrust) or X = z (Drell-Yan) and the brackets denote averaging with 
the corresponding distribution in T or z. For the Drell-Yan cross section it is given by 
eq. ( |3.7| ). The omitted terms do not contribute to terms of order X/Q, since each power 
of 1 — X suppresses the region of phase space X — > 1 where linear terms can arise. One 
can show that (1) is free from such terms in both cases, so that the problem reduces to 
the average of 1 — X. This average receives no virtual corrections and is free from soft 
divergences. 
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For the Drell-Yan cross section, we use the new variable y = 1 — z — £z = X 2 /Q 2 ). 
Neglecting all terms that can not give rise to y/£ in the expansion for small £, we find 



(1 -z) ~ 4 



CfQs 



7T 



dy In 



V + vV - I A Cfois 



TT 



eln£ + 0- >/!+••■ 



(B.2) 



where e is an arbitrary cutoff that does not affect linear terms in A and '~' means that 
the expression is accurate to linear terms. Although the lower limit of the integral is 
its expansion does not contain such terms. The logarithm of £ is due to the fact that the 
partonic cross section is not collinear safe. It is subtracted by the parton distribution 
function, which does not introduce linear terms either (Sect. 3). 

In case of thrust it is convenient to switch from the energy fractions Xi and X2 of 
quark and anti-quark as integration variables to y — 1 — 1/2 (X1 + X2) and x = [x\ — x-i)j2. 
We then find after integration over x 



(l-T) 



CfOLs 



TT 



dy 



CfOL s 



TT 



i-4 



m , = + m 



y - \Jy 2 - i 



O-lnf -4G\/f + 



(B.3) 



where G = 0.91596 ... is Catalan's constant. The additional logarithm in the integrand 
as compared to the Drell-Yan case eliminates the logarithm of £ as required by collinear 
safety of thrust, but introduces a (n additional) linear correction in X/Q. If the in- 
dependence of the normalization J2i \Pi\ of thrust is neglected, the square root in front 
of the logarithm is absent and we obtain the different result of [0. 
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